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Abstract. - The unusually narrow features in the fluorescence from So Rb driven by the 
\Q , cooling and repumper laser fields, reported in [1], are explained on the basis of a four-level 

density matrix calculation. Quantum effects alter the efficiency of atom transfer by the probe 
(repumper) laser to the levels connected by the pump (cooling) laser. This, combined with 
the double resonance condition [1] , leads to velocity selection from co-propagating and counter 
propagating pump and probe beams resulting in narrow fluorescence peaks from a thermal gas 
at room temperature. 

o 

Introduction. - Multilevel atoms, under the action of multiple fields, display a variety 
of phenomena brought about by the ac Stark splitting of electronic levels and by interference 
due to induced coherences. These manifest themselves as splitting of peaks, inhibition of 
absorption or emission, narrow windows of transparency, etc. Theoretical work on this front 
have considered A, V and ladder type of three level systems [2]. The presence of additional 
levels increases multifold the possibilities of interference phenomena that exist. The nested 
A and the N systems have been examined to some extent [3]. In this paper we examine an 
£j \ "inverted N" system using the density matrix formalism. Recent experimental observations 
in S5 Rb [1] have prompted the study of this four-level system under the action of two fields. 

Experiment. The experiment in [1] consisted of driving a collection of 85 Rb atoms 
with two fields- a strong "cooling" field and a weak "repumper" field. One of the lasers, the 
cooling (pump) laser was held at a fixed detuning S c from the conventional cooling transition 
( 1 ) (5Si/2,F — 3 — > 5Ps/2,F' — 4') and the repumper (probe) laser was scanned with its 
detuning (<5 r ) from 551/2, F = 2 — > 5P3/2, F' = 3' level varying over the entire manifold, F=2 
— » F'=l',2',3' (See Fig. 1) ( 2 ). Narrow fluorescence peaks were observed at definite S r for 
a given S c . While one would expect a broad fluorescence, with the width indicative of the 
temperature of the gas, the experiment showed fluorescence peaks just 30MHz in width, much 
smaller than their Doppler width of 500MHz. This was explained on the basis of a double 
resonance model which maximised fluorescence whenever the atom found both the cooling and 
repumping lasers on resonance with their respective transitions. A simple analytical treatment 
based on the above model showed that atom transfer efficiency by the repumper laser to the 

( 1 ) The cooling and repumper beams are used in the sense as defined in [1]. 

( 2 ) Hyperfine levels of 5Si/2 are denoted unprimed while that of 5P3/2 are denoted primed. 
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Fig. 1 - (a): Energy level diagram of 85 Rb. (b):A one dimensional configuration (along Z) of the 
cooling and repumper beams with the S5 Rb atom taken move along +z direction 



levels connected by the cooling laser was maximised under the double resonance condition 
and gave rise to narrow fluorescence peaks. 

In this paper, we present a rigorous treatment of the same problem using a four level 
density matrix formalism. Our treatment differs from the earlier studies of four-level density 
matrix [4] in that (a) it considers an inverted N system and (b) it takes into account the 
motion of the atoms in co- and counter-propagating geometries with respect to the cooling 
and repumper beams. 

Four Level Density Matrix. - The four levels under consideration are: the two ground 
hyperfinc levels F=2,3 and the two excited levels F'=3',4'. For simplicity, we have considered 
an one dimensional situation where the two driving fields are in the ±z direction and the 
atom is moving along the z direction with a velocity v. The detuning of the cooling laser 
taking the Doppler effect into account is A c = 5 C — k c .v and the repumper laser's detuning is 
A r = 6 r — k r .v where, S c and S r are the detunings of the laser in the laboratory frame. The 
total Hamiltonian for the system consisting of the atom and the light fields is written in the 
interaction picture as 

H = H + Hj (1) 

where Hq is the Hamiltonian for the bare atom and Hi is the atom-light interaction Hamil- 
tonian. They are given as 



H = hu 2 \2 >< 21 + fwj 3 \3 >< 31 + hu 3 >\3' >< 3'| + huj 4 ,\4' >< 4'| 



and 



Hi = -~[fi34'|3 >< 4'\exp(-iuj L ct) + fl 33 ,\3 >< i'\exp{-iw L ct) + ^23' exp (-iu LR t)\2 >< i'\+H.C] 

Here the huii represent the energies of the levels as represented in Figure 1, with huj 2 taken 
as zero, lu^c and u>lr the frequencies of the cooling and repumper laser beams and fi^-/ is 
the Rabi frequency connecting the levels i and j' . The total Hamiltonian can be written in 
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where the rows and columns correspond to levels 2, 3, 3', 4' in sequence. The dynamics of the 
system described by this Hamiltonian can be studied using the density matrix p = J2 Pij\i >< 
j\. The time evolution of the density matrix p is given by the Liouville equation 



(2) 



with 



(3) 



where 7i'->j being the spontaneous decay rate from the j level to the i level. ( J ) The 
rate equations of the four levels for an atom moving with a velocity v are derived under the 
rotating wave approximation. They are 



dp 



li 



dt 



"^[^'PSl - ^23'Pl3] + 273'2P33' (4) 



= -«(Ac3' - A r ) - ^[fi 2 3'P32 - ^ 33 'Pl3 ~ ^34'Pl4] (5) 
= (i(A r ) - (73/2 + 73'3))P13 - |[^23'P33 - ^23'PH ~ ^33^12] (6) 



dp 



14 



dt 



= -i[(Ac3/ - A c - A r ) - 74' 3 ]pi4 - 2 [^23'P34' - 3 4'Pi 2 ] (7) 

= i(A c3 / - A r )p 2 l - ^[^33'P31 + ^34'P41 - ^23'P23] (8) 
-^jT~ = _ ^[^33'P32 + ^34'P42 - ^33'/?23 - ^34'P24] + 273*3/333 + 274/3/344) (9) 



(«A c3 / - (73/2 + 73'3))P23 - ^[^33'P33 + ^34'P43 - ^23'/?21 - ^33'P22] 



dt 



dp. 



24 



dt 



(lA c - 74'3)/?24 - 2 ^33'P34 + ^34'P44 - ^34'P22] 



dp 



'31 



dt 

dpZ2 

dt 



= (-iA r - (73/3 + 73'3))P31 - 2 [ fi 23'Pll + ^33'P21 ~ ^23'P3s] 

1 

(-iA c y - (73/2 + 73'3))P32 - 2^33'P22 - ^33^33 - ^34'P3i] 



dp 



'33 



-2(73/2 + 73'3)P33 - 2 [ fi 23'Pl3 + ^33'P23 - ^23'P3l] 



Cfo34 



[(zA c3 ' - A c ) - (73-2 + 73'3 + 74'3)]P34 - ;:[^23'Pl4 + ^33'P24 - ^34'P32] 



dp 41 

dt 
dp. 



42 



dt 



[(iAc3- - A c - A r ) - 74'3]P41 - ^ [^34P21 - ^23'P43] 
[(iA c - 74' 3 )]P42 - 2^34^22 ~ ^ 33 'P43 - ^34'P44] 



(10) 

(11) 
(12) 

(13) 
(14) 
(15) 
(16) 
(17) 



( 3 )where F^i is the decay rate of the level j to the level i'. This is however a diagonal matrix 
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Fig. 2 - Trace 1 shows the saturation absorption spectrum of the repumper. Trace 2 shows the 
experimental curve (black) and the corresponding theoretical curve (blue). See text for explanation 
on the (a) and (b) peaks in the figure. 

;[^34'P23 - ^33'P42 ~ ^3 4 'P23'P4l] (18) 
-2^3i'P^ - ^34'P42 - 274-3/544] (19) 

For each velocity v steady state values of pij are obtained by numerically solving the above 
rate equations for various values of S c and S r , subject to the constraint Yli=i Pa ~ Thus 
for an atom with a velocity v we obtain the population of each level and coherence between 
various levels for different values of 8 r and S c . The fluorescence emitted by atoms with a 
velocity v is given by 

Fluorescence(A c , A r ) = r 4 /3p4/4/ + r 3 / 3 p33/ + r 3 / 2 p3'3' (20) 

As the detector collects fluorescence from atoms which are in thermal motion, we take the aver- 
age of each p^ value over the range of velocities, weighted by the one dimensional Maxwellian 
velocity distribution. The fluorescence calculated as above as a function of S r and the corre- 
sponding experimental data are given in Fig. 2 for a cooling laser detuning S c — —162 MHz. 
A general agreement between the results of our calculation and that of the experiment is seen. 
( 4 ) The individual features will be discussed in detail below. 

Discussion. - Consider a situation shown in Fig. lb, when the cooling and repumping 
beams are both along ±z directions. Initially let us consider the atoms to be at rest. For a 
given detuning S c of the cooling (pump) laser, we should get fluorescence peaks corresponding 
to the Autler-Townes (AT) dressed states of F' = 3' at the probe (repumper) detunings [5] 

S r± = 6 -f±^/(6i 3 ,+W)/2 (21) 

( 4 ) Peaks (a) and (b) seen in the experiment arise due to the transitions from F 1 = 2' level which has not 
been included in our density matrix 



dp43 
dt 



= H(A c3 ' - A c 



(73' 2 + 73'3 + 74'3)]P43 " 

dt 
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Fig. 3 - The imaginary part of P23' vs the positions of the probe absorption S r showing that absorption 
of the repumper (probe) light due to Autler -Townes splitting of the level F' — 3' occurs at the 
detunings 0MHz and — AlMHz. (b) P33 (blue) and P22 (green) as functions of 5 r - 

Here S c y = 2ir* 121MHz + 8 C denotes the detuning of the cooling laser from 3^3' transition, 
121 MHz being the level spacing between F' = 3' and 4' levels and f2 its Rabi frequency. For 
the specific case shown in figure 2, for a cooling laser detuning S c = —l&2MHz and for small 
fl, the Autler- Townes peak positions for these zero velocity atoms are S r ± rj 0, —41 MHz. 
This is confirmed from our density matrix calculations. At these detunings of the probe its 
absorption is maximised and repumping most efficient as seen from Figs. 3(a) and 3(b) which 
give P23', P33 and P22 as functions of 8 r . 

Consider now an atom in motion with a velocity +v, along the +z direction. If S c < 
this atom predominantly absorbs from the cooling beam coming towards it (C). Since the 
repumper laser is scanned, depending on the sign of S r the atom absorbs either from B or 
from D. So for a given 8 C , altogether we get four peaks, a Autler- Townes pair each shifted 
due to the Doppler effect, depending upon whether absorption takes place from the repumper 
beam B or from repumper beam D. The same holds for S c > for atoms with velocity — v . 

The AT doublet positions as given in equation (21) when applied to atoms in motion and, 
for small f2 and S c < are 

<5 r _ - (±kv) « 0; 8 r+ - (±kv) » S c3 , - (±kv) (22) 

As the repumper is scanned, maximum transfer of population from F =2 to F= 3 will occur 
when these conditions are satisfied. Now, as mentioned in [1] and as will be shown shortly, 
the fluorescence is predominantly from 4'— > 3. However, not all atoms in F =3 will be in 
resonance with the cooling laser; only a small velocity class around v c = 8 c /k will give rise 
to fluorescence. This is the mechanism that gives rise to narrow velocity selection from a hot 
gas. This velocity selection effect is confirmed by our calculations and is shown in Fig. 4 
which gives and pyy as functions of the velocity of the atom for S c = S r = -162 MHz. We 
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Fig. 4 - Populations in the upper levels F' = 3 and F' =4'. The very low population of F' — 3' in 
(a) is shown magnified in (b). 



see that the population in F'=4' (pw) is two orders more than that in F'=3' (p3'3') and it 
peaks at the critical velocity v c = 126m/ sec. Therefore only for velocities around v c will the 
fluorescence from F' = 4' be maximised, showing that indeed the double resonance condition 
selects a narrow velocity class for fluorescence. 

The Autler-Townes peak positions, as given by equations ( 12211 . with the constraint that 
double resonance should also be satisfied for maximum fluorescence, require that 

<5 r _ = S c - 6 r+ = 121MHz + S c (23) 

when repumper absorption takes place from D and 

<5 r _ = -S c ; S r+ = 121MHz - S c (v = v c ); S r+ = 121MHz + S c + 2kv(v ^ v c ) (24) 

when repumper absorption takes place from B. 

The same conditions result for atoms with a velocity —v, for S c > 0. It has been estimated 
in [5] that the peak at S r - is broad and the one at 5 r + is narrow. The width of the AT peaks 
decide the prominence of a fluorescence peak as it decides the number of atoms participating 
in the fluorescence. Thus for peak positions 5 r - = 5 C & —6 C we expect a large fluorescence 
whereas for the peak at 5 r + = 121MHz + 5 C we expect a much smaller fluorescence as is 
indeed seen from experiment and from our calculations (Fig. 2). As mentioned in [1] this 
peak is well resolved only at large detunings. 

The peak around <5 r + = 121 MHz + S c +2kv is absent both in experiment and in our density 
matrix calculation. For this case, the repumper absorption takes place from B whereas the 
cooling is absorbed from C. When the absorption takes place from counter-propagating cooling 
and repumper beams the velocity class satisfying the double resonance is severely restricted. 
In fact only for v = v c will the double resonance condition be satisfied. Atoms with v ^ v c 
will see the cooling and repumper to be shifted by different detunings and hence will not 
contribute to the fluorescence. Thus the peaks resulting from this configuration will not be 
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Fig. 5 - Population in the level F' — 4' as a function of the velocity of atoms for the Autler-Townes 
peaks at (a) 5 r = -5 C = 162MHz (b) 5 r = V2\MHz + 8 C + 2kv. 



resolvable as only a very small number of atoms contribute to it. The severely restricted 
velocity range at resonance resulting in the absence of the peak at \2\MHz + 5 C + 2kv is 
illustrated in Fig. 5b which is obtained from the density matrix calculation. As we see here, 
only a small number in a very narrow velocity range contribute to population in F'=4' in 
contrast to the case 6 r _- = — S c (Fig. 5a). 

The peaks marked (a) and (b) in Figure 2 are the peaks corresponding to the AT levels 
of F' = 2' at <5 r _ = S c & — S c corresponding to absorption from counter and co-propagating 
repumper beams. Our calculation has not reproduced this peak, as the level F'=2' was not 
included in the density matrix. Nevertheless the peak positions can be simply calculated from 
the corresponding peaks for F'=3' by shifting them by —Q3MHz. The peak at 8 C + 121MHz 
falls at the same 5 r + for both F' = 2' and F' = 3'. 

The theoretically calculated widths of the fluorescence peaks match the experimentally 
observed narrow width of about 30MHz. It should be emphasised that this does not arise 
due to the cooling of atoms in the optical molasses like configuration but due to the velocity 
selection as discussed above. These fluorescence peaks are experimentally seen to be narrow 
even for blue detunings of the cooling beam where no cooling occurs and the atoms are at room 
temperature. The corresponding Doppler width is several hundred MHz at this temperature. 
The additonal effects of cooling on these fluorescent peaks in a Doppler free beam geometry 
will be discussed elsewhere. 

Conclusions. The experimental observations of [1] has been explained using a four- 
level density matrix formalism. The theory finds that the fluorescence peaks are given rise 
by the Autler-Townes (AT) doublets of F' = 3' and F' = 2' for atoms around a particular 
velocity class. The theory gives the three AT peaks of F' = 3' at positions borne out by the 
experiment. The theory also explains the absence of the fourth peak. The widths and heights 
of the peaks agree quite well with the experimentally obtained widths and heights. 
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